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ABSTRACT 

Consideration  is  given  to  a  system  of  reaction  diffusion  equations  which 
have  qualitative  significance  for  several  applications  including  nerve 
conduction  and  distributed  chemical/biochemical  systems*  These  equations  are 
of  the  FitzHugh-Nagumo  type  and  contain  three  parameters.  For  certain  ranges 
of  the  parameters  the  system  exhibits  two  stable  singular  points.  A  singular 
perturbation  construction  is  given  to  illustrate  that  there  may  exist  both 
pulse  type  and  transition  type  traveling  waves.  A  complete ,  rigorous, 
description  of  which  of  these  waves  exist  for  a  given  set  of  parameters  and 
how  the  velocities  of  the  waves  vary  with  the  parameters  is  given  for  the  case 
when  the  system  has  a  piecewise  linear  nonlinearity.  Numerical  results  of 
solutions  to  these  equations  are  also  presented.  These  calculations 
illustrate  how  waves  are  generated  from  initial  data,  how  they  interact  during 
collisions,  and  how  they  are  influenced  by  local  disturbances  and  boundary 
conditions.  In  one  example,  a  rightward  moving  "front"  (rest-to-excited 
transition)  slows  down,  stops,  turns  around,  and  develops  into  a  leftward 
moving  "back"  ( excited-to-rest  transition)  when  it  encounters  appropriate 
spatially  localized,  transient  perturbations  in  the  dependent  variables.  In 
this  case,  the  medium  ultimately  returns  to  the  rest  state.  For  the  same 
parameter  values  we  have  sought,  but  have  not  found,  analogous  local 
disturbances  sufficient  to  reflect  a  propagating  "back".  These  observations 
are  consistent  with  the  notion  that  the  rest  state  is  in  some  sense  dominant 
for  a  certain  parameter  range*  It  is  conjectured  that  local  disturbances  or 
certain  boundary  conditions  cannot  cause  the  reflection  of  a  nondominant  to  a 
dominant  transition  wave.  Partial  evidence  for  this  conjecture  is  presented 
by  means  of  an  analytic  argument. 
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SIGNIFICANCE  AND  EXPLANATION 


The  mathematical  equations  studied  here  were  originally  introduced  as  a 
model  for  nerve  conduction*  They  have  also  been  considered  as  a  model  for 
chemical/biochemical  reaction  systems  as  well  as  electrical  transmission 
lines*  It  has  been  demonstrated  that  electrical  impulses  in  the  nerve  axon 
travel  with  constant  shape  and  velocity*  Mathematically  this  corresponds  to  a 
traveling  wave  solution*  For  a  given  mathematical  model  there  may  exist 
different  types  of  traveling  wave  solutions  each  traveling  with  a  different 
velocity*  These  include  solitary  pulses,  multiple  pulses,  and  traveling 
fronts. 

In  this  report  we  describe  which  waves  exist  for  a  given  set  of 
parameters  of  the  equations  and  how  the  velocities  of  the  waves  vary  as 
functions  of  the  parameters*  Numerical  results  of  solutions  to  these 
equations  are  presented  to  illustrate  how  the  waves  are  generated  from  initial 
stimuli,  how  they  interact  during  collisions,  and  how  they  are  influenced  by 
local  disturbances. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


PROPAGATION  PHENOMENA  IN  A  BISTABLE  REACTION  DIFFUSION  SYSTEM 


John  Rinzel  and  David  Terman 


1.  Introduction* 

The  FitzHugh-Nagumo  (FHN)  equations 

v  -  v  +  f(v)  -  w 
t  XX 

(1.1) 

wt  *  e*v  "  ^  '  e,Y  >  0 

has  been  studied  extensively  as  a  qualitative  model  for  nerve  conduction  [7,  12,  18).  In 
this  context  the  dynamics,  (1.1)  with  vxx  =  0  ,  are  usually  taken  to  be  excitable  with 
nullclines  as  in  Figure  1A.  The  rest  state  (0,0)  is  stable  and  globally  attracting  but 
an  adequate  initial  displacement  from  rest  leads  to  a  large  excursion  before  the  eventual 
return  to  rest.  Qualitatively  similar  dynamics  are  found  for  other  excitable  systems  and 
study  of  the  FHN  equation  has  provided  insight  into  excitability  and  propagation  phenomena 
in  several  applications,  e.g.  see  [18). 

Here,  we  extend  consideration  to  parameter  ranges,  e.g.  y  large  enough,  for  which 
the  v  -  w  dynamics  are  bistable  (Figure  IB).  In  this  case  there  are  two  stable  singular 
points:  the  rest  state  on  the  left  branch  of  f(v)  and  the  excited  state  (EV,EW)  on  the 
right  branch.  Although  this  parameter  range  precludes  direct  applicability  to  most  nerve 
or  muscle  membrane  types  there  are  examples  of  conditions  which  induce  bistability  (8, 

11] •  Dynamics  similar  to  Figure  IB  have  also  been  considered  theoretically  in  models  for 
neuronal  interactions  at  the  population  level  [19].  For  chemical/biochemical  reaction 
systems  [14],  as  well  as  electronic  transmission  lines  [13),  two  variable  systems  with 
bistable  behavior  have  been  investigated.  Sometimes  in  these  studies,  although  not 
necessarily,  the  w-variable  has  much  slower  dynamics  than  v  in  which  case  0  <  e  <<  1. 
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This  spatial  coupling  is  physically  natural  for 


2 

In  (1.1),  3^  operates  only  on  v  . 

the  electronic  transmission  line  as  well  as  the  nerve  fiber.  For  the  neuronal  population 
models,  this  corresponds  to  short-range  coupll-ng  of  one  cell  type  (say,  the  inhibitory 
cells)  relative  to  the  length  scale  of  coupling  for  the  other  cell  type  (the  excitatory 
cells).  For  chemical  systems,  one  might  neglect  wxx  if  w  ,  the  species  without  the 
cubic-like  instantaneous  dynamics,  diffuses  much  more  slowly  than  the  v-species.  In  (14], 
for  an  enzyme- media ted  reaction  scheme,  the  assumption  of  very  rapid  v-dynamics  is  used  to 
scale  the  model  equations  and  obtain  a  system  qualitatively  like  (1.1)  and  Figure  IB  with 
only  v  as  diffusible. 

In  this  paper  we  will  investigate  various  propagation  phenomena  for  (1.1).  An  example 
of  Ortoleva  (described  by  Fife  (6])  offers  preliminary  motivation.  Suppose  that  at  t  *  0, 
the  medium  is  in  the  excited  state  for  x  >  0  and  in  the  rest  state  for  x  <  0.  Subse¬ 
quently,  a  transition  wave  will  develop  which  then  propagates  away  from  x  *  0  so  that  the 
entire  medium  is  ultimately  brought  to  one  of  the  two  stable  states,  whichever  is  dominant 
(a  descriptive  term  also  used  in  (6,  14]).  Hence  (1.1)  should  have  a  traveling  wave 
solution  which  corresponds  to  a  steadily  propagating  transition  either  from  the  rest  state 
to  the  excited  state,  or  vice  versa,  (Figure  2  A,B)  depending  on  which  state  dominates.  As 
it  turns  out,  for  a  certain  parameter  range,  both  types  of  transition  waves  exist,  simulta¬ 
neously  with  different  speeds;  we  refer  to  these  as  traveling  front  and  back  solutions. 
Intriguingly,  for  a  subset  of  the  front-back  parameter  range,  traveling  pulse  solutions  are 
also  found.  A  pulse  represents  an  excursion  either  from  the  rest  state  back  to  itself  or 
from  the  excited  state  back  to  itself  (Figure  2  C,D).  In  addition,  there  may  be  solutions 
for  steadily  propagating  trains  of  pulses  although  we  will  not  treat  trains  here. 

In  Section  2  we  will  outline  the  singular  perturbation  construction  [2,  3,  91  of  the 
various  traveling  waves  for  0  <  e  <<  1  .  This  will  develop  our  intuition  for  the  parameter 
ranges  in  which  the  waves  exist.  Two  special  cases  of  (1.1)  are  treated  explicitly:  CUB, 
when  f(v)  =  v(v-a)(1-v)  and  PWL,  when  f(v)  -  -v  +  H(v-a)  where  H( *)  is  a  Heaviside 
step  function?  here,  0  <  a  <  V2  *  Next,  in  Section  3,  we  obtain  analytically  the  front, 

back,  and  pulse  solutions  for  PWL  in  a  large  parameter  range  which  extends  beyond  the 
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reqion  0  <  e  <<  1.  Rigorous  proofs  for  existence  of  these  waves  are  given.  We  also 
calculate  wave  speeds  as  functions  of  the  parameters.  In  addition  to  stable  waves,  we  find 


unstable  ones.  A  global  picture  (Figure  6)  is  provided  for  e, ^-parameter  regimes  in 
%  which  the  various  waves  exist. 

Finally,  in  Section  4,  we  present  numerical  results  of  solutions  to  (1.1),  PWL,  for 
given  initial  and  boundary  conditions.  These  calculations  illustrate  how  waves  are 
generated  from  initial  data,  how  they  interact  during  collisions,  and  how  they  are 
influenced  by  local  disturbances.  In  one  example,  a  rightward  moving  front  (rest-to- 
excited  transition)  slows  down,  stops,  turns  around,  and  develops  into  a  leftward  moving 
back  ( excited-to-rest-transition)  when  it  encounters  appropriate  spatially  localized, 
transient  perturbations  in  v  and  w  «  In  this  case,  the  medium  ultimately  returns  to  the 
rest  state.  For  the  same  parameter  values  however,  we  have  sought,  but  have  not  found, 
analogous  local  disturbances  sufficient  to  reflect  a  propagating  back.  These  observations 
are  consistent  with  the  notion  that  the  rest  state  is  in  some  sense  dominant  for  these 
parameter  values.  We  conjecture  that,  for  a  certain  parameter  range,  local  disturbances 
may  lead  to  reflection  of  a  front  but  cannot  cause  reflection  of  a  back,  or  more  generally, 
cannot  cause  reflection  of  a  non-dominant  to  a  dominant  transition  wave.  We  offer  partial 
support  for  this  conjecture  with  an  analytic  argument.  Similarly  we  expect  for  such 
parameters  that,  when  a  boundary  is  held  at  rest,  an  approaching  front  will  be  reflected. 

On  the  other  hand,  a  back  which  propagates  toward  a  boundary  held  at  the  excited  state  will 
not  reflect  but  rather  may  tend  to  a  spatially  nonuniform,  t- independent ,  steady  state,  a 
solution  which  decays  from  the  excited  to  the  rest  state  with  distance  from  the  boundary. 
For  a  chemical  system,  these  phenomena  would  be  observed,  for  example,  in  a  capillary  tube 
which  makes  contact  at  the  boundary  with  an  infinite  bath?  where  the  bath  maintains  either 
the  rest  or  excited  state. 

The  influence  of  boundaries  and  initial  disturbances  on  propagating  transitions,  and  a 
notion  of  dominance,  have  also  been  discussed  by  Fife  [6].  In  his  (as  well  as  others) 
singular  perturbation  treatment,  fronts  are  formed  because  of  pseudo-steady  state 
bistability;  that  is,  the  well-studied  reduced  equation,  (1.1)  with  e  -  0  and  w  , 
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constant,  axhibits  a  transition  wavs*  In  such  cases,  only  one  type  of  stable  transition 
wave  exists  and  the  definition  of  dominance  follows  naturally;  the  wave  is  always 
monotone.  For  the  full  system  (1.1),  fronts  and  backs  are  generally  not  monotone  but 
rather  have  over  or  undershoots  (Figure  2  A,B)  or,  in  some  cases,  damped  oscillatory 
approach  to  the  steady  state  (also,  see  [15)).  The  bistability  and  (not  yet  rigorously 
defined)  notion  of  dominance,  which  we  discuss,  is  for  the  full  system.  Fronts,  backs,  and 
pulses,  were  also  considered  by  Collins  and  Boss  [4],  Keener  [9],  and  Hagumo,  et  al  (13); 
Ortoleva  and  Ross  [14)  describe  fronts  and  backs,  but  not  pulses,  for  a  bistable  system. 
KLaasen  and  Troy  [10]  also  treat  some  behavior  of  two-variable,  bistable  systems. 


2.  The  Singular  Perturbation  Construction  of  Traveling  Waves* 

In  this  section  we  describe  ways  of  determining  intuitively  which  traveling  waves 

exist  for  a  given  set  of  parameters  and  how  their  velocities  change  as  functions  of 

parameters.  Our  description  summarizes  the  results  of  (2], [3],  [9] •  If  { vQ( z) , wc ( z) ) , 

z  *  x  +  ct,  is  a  traveling  wave  solution  of  system  (1.1),  then  setting  u_(z)  *  v'(z)  , 

c  c 

it  follows  that  Vc(z)  -  (vc(z),  uc(z),wc(z))  satisfies  the  first  order  system  of  ordinary 
differential  equations 


v'  =  u 
c  c 


(2,1) 


cu  -  f(v  )  +  w 
c  c  c 


( V 


yw  ) 


c  c  c 

The  number  of  equilibrium  states  of  system  (2,1)  will  be  either  one,  two,  or  three 

depending  upon  how  the  line  v  =  yw  intersects  the  cubic-like  curve  w  =  f(v  ), 

c  c  c  c  ' 

u  =0.  if  the  parameters  a  and  e  are  fixed,  then  when  y  is  sufficiently  small  the 
c 

origin  will  be  the  only  equilibrium  state  of  system  (2.1).  (Figure  1A).  For  y  large 

there  will  be  three  equilibrium  states.  In  this  case  we  shall  refer  to  the  origin,  as 

the  rest  3tate  and  E  *  (Ev,0,Ew)  as  the  excited  state  (Figure  IB).  We  shall  assume 

c  >  0  so  Vc  corresponds  to  a  leftward  moving  wave;  note  that  the  spatial  profiles  in 

Figure  2  are  for  the  rightward  moving  versions. 

For  a  pulse  we  seek  a  value  of  c  such  that  system  (2.1)  has  a  solution  which 

satisfies  lim  V  (z)  *  0,  a  front  is  a  solution  of  (2.1)  which  satisfies  iim  V  (z>  = 

|  z  |  -h»  z+-°° 

and  lim  V  (z)  =  E,  while  a  back  satisfies  lim  V  (z)  *  E  and  lim  V  (z)  =  0,  In 

Z>+<»  J5+—  oo  £++« 

the  case  of  multiple  equilibrium  states  it  is  also  possible  to  have  solutions  which  satisfy 
lim  V(z)  *  E.  We  shall  refer  to  these  as  E -pulses . 


(z  |  >oo 

We  now  sketch  the  singular  perturbation  construction  of  traveling  waves.  The  solution 
of  system  (2.1)  is  sought  in  separate  regions,  labelled  "inner"  and  "outer",  and  then 
joined  appropriately.  The  inner  region  is  defined  to  be  those  intervals  over  which  rapid 
changes  in  vc  take  place,  while  the  outer  region  involves  variations  on  a  much  larger 
z  scale. 


0 
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We  first  make  the  change  of  variables  £  *  ez  in  (2.1)  to  obtain: 


eu'  -  cu  -  f(v  )  +  w 
c  c  c  c 


(  *»  ‘at 


w*  *  '  (v  -  yw  )  ; 
c  c  c  c 


Outer  Region.  Away  from  the  locations  of  sharp  transitions  in  vc  we  obtain  the  lowest 
order  approximation  by  setting  e  *  0  in  (2.2)  to  obtain: 


w  »  f(v  ) 
c  c 


(  #)  "  ‘ 


-  (v  -  yw  ) 
c  c  c 


This  implies  that  in  the  outer  region  the  solution  lies  on  the  slow  manifold  curve: 

wc  *  f<vc),  u^  =  0  .  We  assume  that  the  left  branch  of  the  cubic  i9  given  by  v  *  g"(w) 


and  the  right  branch  is  given  by  v  -  g*(w). 


Inner  Region.  Here  we  use  the  stretched  variable  z  -  — £  and  system  (2.2)  returns  to 

system  (2.1).  With  e  *  0  in  (2.1)  we  obtain  the  following  lowest  order  approximation 

f  v*  *  u 
c  c 


u*  “  cu  -  f(v  )  +  w  , 
c  c  c  c 


<•)  -  —  . 
dz 


Therefore,  in  the  inner  region  the  solution  must  lie  in  a  plane  w  =  constant.  We  shall 

c 

see  that  this  constant  determines  the  speed  of  the  traveling  wave. 

Matched  Composite  Solutions.  To  lowest  order  in  e  »  the  full  solution  is  obtained  by 
matching  the  inner  and  outer  solutions  where  their  domains  of  definition  overlap.  We  treat 
each  wave  separately. 

Fronts.  We  assume  that  the  parameters  are  chosen  so  that  the  line  v^  *  Ywc  intersects  the 
cubic  as  shown  in  Figure  IB.  Since  the  independent  variable  £  does  not  appear  explicitly 
in  Equation  (2.2)  we  may  also  assume  that  the  inner  region  lies  in  a  small  neighborhood 


about  the  origin,  £  *  0 


The  outer  region  describes  the  front's  plateau.  For  Equation  (2.3)  we  take  vc  =  g+(wc) 
and  the  initial  condition: 

(2.5)  wc(0)  =  0  . 

This  means  that  w  *  E  as  z  «. 

c  w 

In  order  to  match  the  inner  and  outer  solutions  we  take  the  following  boundary 
conditions  in  Equation  (2.4): 

(2.6)  v  (-«)  -  0  ,  v  (  +  «)  *  1  ,  w  (z)  =  0  . 

c  c  c 

Equation  (2.4)  for  cubic-like  f(v)  with  boundary  conditions  (2.6)  was  studied  by  Aronson 

and  Weinberger  (11*  They  proved  the  existence  of  a  unique  solution  for  a  unique  value  of 

c  provided  A  >  0  ,  where  A  =  /^ffvjdv.  Hence,  these  equations  determine,  up  to  lowest 

0 

order  in  e  ,  the  speed  of  the  front.  In  the  case  f(v)  given  by  (PWL)  one  can  give  an 
explicit  formula  for  the  speed,  c,  of  solutions  of  Equation  (2.4)  with  boundary  condition 

(2.6)  as  a  function  of  the  parameter  a  (see  [16]).  This  formula  is 

c  -  (1  -  2a) [a(1  -  a)]"  ^2  . 

Note  that  c  is  a  decreasing  function  of  a  and  lim  c  *  0  .  This  is  similarly  true 

at  y2 

for  f(v)  given  by  (CUB),  in  which  case  c  *  /2  (  V2  -  a)  (see  [3]). 

Note  that  the  outer  solution  for  a  front  lies  on  the  right  branch  of  the  cubic  while 
the  inner  solution  lies  in  the  plane  w^  =  0  .  The  complete  phase  space  trajectory  for  a 
lowest  order  front  is  shown  in  Figure  3. 

Backs.  In  order  for  a  back  to  exist  the  line  v  -  yw  must  intersect  the  cubic  as  shown 
-  c  c 

in  Figure  IB.  Again  we  may  assume  that  the  inner  region  lies  in  a  small  neighborhood  about 
K  -  0*  For  the  outer  region,  we  take  vc  =  g~(w  )  in  Equation  (2.3)  and  the  initial 
condition : 

(2.7)  wc(0>=£w- 

In  this  case,  w^.  <  0  and  w  -►  0  as  z  ♦  ».  Then  in  Equation  (2  4)  we  take  the  boundary 
conditions 

(2.9)  v  (-«)  =  E  ,  v  (+»)  =  g~ ( E  )  ,  w  (z)  H  5  . 

c  V  c  w  c  w 


The  results  of  Aronson  and  Weinberger  imply  that  there  exists  a  unique  solution  of  Equation 


«n  the  right  branch  to  the  left  branch  of  the 
he  trajectory  of  a  pulse  occurs  in  the  plane 

so  that  the  shaded  regions  have  equal  areas. 


(2.4)  with  boundary  conditions  (2.8)  for  a  unique  positive  value  of  c  if  and  only  if 


,« (EW> 

B  =  {  [f(v)-E]dv>0.  From  this  we  conclude  that  the  back  (to  lower  order)  does 

*  w 

£v 

not  exist  if  y  is  too  large.  This  is  because  the  assumption  A  >  0  implies  that  there 
must  exist  a  critical  value  of  y  ,  y  =  ,  such  that  B  <  0  for  y  >  y2  .  Hence,  the 

"singular"  or  lowest  order  back  cannot  exist  for  y  >  y^  •  For  f(v)  given  by  (PWL)  a 

1+2a 

simple  computation  shows  that  y2  =  -  *  In  this  case  one  can  also  give  an  explicit 

formula  for  the  speed,  c,  of  solutions  of  Equation  (2.4)  with  boundary  conditions  (2.8) 
as  a  function  of  y  and  a  .  This  formula  is 

c  •  t2a  "  ><a  '  TT7  1 1/2  • 

We  see  that  for  the  parameter  a  fixed  the  speed  is  a  decreasing  function  of  y  ,  and 
lim  c  «  0  .  These  properties  also  hold  for  f(v)  given  by  (CUB). 


Note  that  the  outer  solution  of  a  back  lies  on  the  left  branch  of  the  cubic-like  curve 
while  the  inner  solution  lies  in  the  plane  w^  =  .  The  total  phase  space  tajectory  for 

a  back  is  shown  in  Figure  3. 

Pulses .  A  pulse  can  be  thought  of  as  a  wave  resembling  a  front  folLowed  by  a  back 

traveling  at  the  same  speed.  The  trajectory  in  the  phase  space  for  the  front  part  of  the 

pulse,  or  upstroke,  begins  at  the  origin  and  then  makes  a  fast  transition  in  the  plane 

w  =  0  to  the  right  branch  of  w_  “  f(v_)  .  It  then  slowly  moves  up  the  right  branch  of 
c  c  c 

the  cubic-like  curve  until  it  reaches  some  point  (Vq,0,Wq);  this  forms  the  plateau  of  the 
pulse.  The  trajectory  for  the  back  part  of  the  pulse,  or  downstroke,  begins  at  the  point 
(v0,0,Wq),  makes  a  fast  transition  in  the  plane  w  =  wQ  to  the  left  branch  of  wc  =  f(vc), 
and  then  slowly  returns  to  the  origin  along  this  left  branch  to  form  the  recovery  phase. 

From  Equation  (2.4)  it  follows  that  the  velocity  of  the  back  part  of  the  pulse  is 
determined  by  the  constant  wQ  .  We  must  choose  wQ  so  that  the  velocities  of  the  front 
and  back  parts  of  the  pulse  are  equal.  If  f(v)  is  given  by  (CUB)  or  (PWL)  then  our 
previous  discussions  imply  that  must  be  chosen  so  that 
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For  f(v)  given  by  (CUB),  choose  yQ  so  that  the  line  v  »  yQw  intersects  the  local 

maximum  of  the  cubic  as  shown  in  Figure  4.  Choose  y^  and  so  that,  respectively, 

the  horizontally  and  vertically  shaded  regions  in  Figure  4  have  equal  areas.  For  f(v) 

given  by  (PWL)  choose  y  =  ,  y  *  7^-  ,  and  y  =  ?  +  .  The  singular  perturbation 

0  1 —  a  1  i— a  2  1 —  2a 

construction  then  suggests  that  the  following  proposition  is  true. 


Proposition  2.1.  For  0  <  e  <<  1  a  pulse  exists  for  y  e  (0,y^)  ,  a  front  exists  for 

ye  (yQ,«>)  •  and  a  back  exists  for  ye  • 

Note  that  all  three  types  of  waves  exist  for  ye  (y^y^)*  We  saw  in  the  singular 

perturbation  construction  that  the  speed  of  the  front  is  ,  to  lowest  order  in  e  ,  a 

decreasing  function  of  A  =  J^f(v)dv  ,  while  the  speed  of  the  back  is  a  decreasing 

0 

9_(E  ) 

function  of  B  »  f  If(v)  -  E  )dv  .  Note  that  if  ye  ( y_  ,y, )  then  A  >  B  .  Hence, 

J  w  0  1 

v 

when  all  three  types  of  waves  exist,  the  back  is  faster  than  the  front.  For  ye  (y<J,y2>, 
the  front  is  faster  than  the  back,  and  in  this  parameter  range  an  E-pulse  exists. 
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It  Is  of  course  possible  for  other  waves  to  exist  which  do  not  arise  from  the  singular 
perturbation  construction.  We  shall  see  in  the  next  chapter  that  there  may  exist  two 
traveling  waves  of  a  particular  type  for  a  given  set  of  parameters.  The  singular 
perturbation  construction  gives  the  faster  wave  which  one  expects,  from  numerical  evidence, 
is  stable.  The  slower  wave  is  expected  to  be  unstable. 
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Solutions  of  the  Piecewise  Linear  Model 


In  this  section  we  consider  system  (1.1)  with  f(v)  defined  by  (PWL)  and  we  give  a 
complete,  rigorous  description  of  which  traveling  waves  exist  for  a  given  set  of  parameters 
and  how  their  speeds  vary  with  the  parameters.  The  construction  we  use  here  follows  the 
work  of  Rinzel  and  Keller  (16). 

Pulses:  We  seek  values  of  c  >  0  for  which  system  (2.1)  has  a  pulse  shaped  solution. 
Because  the  independent  variable,  z  ,  does  not  appear  in  the  equation  explicitly  we  can 
choose  the  origin,  z  *  0  ,  so  that  vc<0)  *  a  .  We  also  assume  that  vc(z^)  *  a  for 
some  Zy  >  0.  Therefore,  the  solution  we  seek  is  a  leftward  moving  version  of  the  form 

illustrated  in  Figure  2 C,  to  be  distinguished  from  an  E: -pulse  (Figure  2D).  It  follows  that 

T 

along  with  Equation  (1.1)  the  solution  Vc  *  (vc,uc,wc)  must  satisfy  the  conditions 


v«  r  -  -i 

c  o 


vc(°)  *  vc(z1 )  =  a  '  z1  >  0 


This  jump  condition  (3.1)  results  from  the  discontinuity  of  f(v).  To  represent  a  pulse 
the  solution  must  also  satisfy  the  condition  v^(z)  +  0  as  I2 I  *  00  • 

We  express  the  solution  as 


3  °SrZ 

U*kV  +v 


I  ck  \  e 

k-1  K  K 


0  <  z  <  z. 


Zy  <  Z 


where  ,  c^,  are  the  roots  of  the  characteristic  polynomial 


*3  +  (£H£-)X2  .  UY+1)X  .  (CCS)  , 


and  where 


. . mar-  I  Vit  *v  I  Vi T if  rfr  r*Tfr  V'  mfrer* 


-JL 

1+y 

1 

0 

and  “ 

1 

-1 

k  1+Y  , 

e(c<x  ♦  ey) 

k  -  1,2/3* 


From  (3.4)  it  follows  that  either 


°1>0,  a3<ol2<0 


«1  >  0  ,  ,  Rea2 


Because  v  +0  as  z  *  •  we  must  take  c«  •  0  •  It  is  this  last  condition  which  will 
c  < 

give  us  a  relationship  for  the  speed,  c,  in  terms  of  the  parameters  a,e,  and  y. 

Using  the  continuity  and  jump  conditions  at  z  ■  0  we  obtain 


„  6l  «  _{2  «  '53 

B-  ■  a  -  — r  ,  B*  ■  — : —  ,  B  rn  —J— 

1  p;  2  3  fa 


Matching  at  z  *  z^  yields 


eir  * 

^  ’  oak 


and  -  p'  (  0^) 


Vi  ^  5i  M  Vi 

c  •  ae  +  —  ( 1  -  e  ) 
P1 


Since  c^  ■  0  ,  we  obtain 


'Vi  M  pi 

e  »(1-a-=->5s 


Using  (3.8)  and  the  fact  that  vc(Zj)  *  a  we  are  led  to  the  following  transcendental 
equation: 


-m* f 
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(3.9) 


F(a, e,c,y) 


=  2-s  +  v"  —7  s 
51  P2 


52  p;  _a2/ai ,  s p;  -“3 


/a. 


61  P3 


o264 


0  . 


2  1  +  Y 

Here  we  set  o  =  — ^-L  •  Note  that  as  a  function  of  s  Equation  (3.9)  does  not  depend  on 


the  parameter  a  •  Equation  (3.8)  gives  a  in  terms  of  s  . 

The  following  theorem  is  proved  in  Appendix  A  by  analyzing  the  roots  of  Equation 
(3.9). 

Theorem  3.1:  Assume  that  the  parameters  e  and  y  are  given. 


(i) 


There  is  a  value  of  the  parameter  a  for  which  a  pulse  exists  if  and  only  if 


ty 2  <  1 


(ii) 


For  a  given  value  of  c  at  most  one  pulse  can  exist  for  a  uniquely  determined 


value  of  a  • 
2 


(iii)  If  ty  <  1  then  positive  constants  cm±nt  cmax  exist  with  cm^n  <  cmax  and  such 
that 


(a)  for  each  c  e  ^cmin#cmax^  there  is  a  unique  value  of  a  for  which  a 
pulse  exists. 


(b)  lim  a  3  0  , 


etc 


min 


(c)  lim  *  0  where  z1  was  defined  by  v(z^)  =  a,  Zy  >  0  , 


ct -c 


min 

(d)  lim  a 


etc 


2 ( 1+y)  ' 


max 


(e)  lim  z. 


etc 


(f)  as  ty  t  1  both  c  , ^  and  c _  t  0  . 

min  max 


In  the  proof  of  this  result  we  will  also  see  that  a  pulse  cannot  exist  if  c  is 
either  too  small  or  too  large. 

We  have  solved  Equation  (3.9)  numerically  for  fixed  values  of  t  ,  c  and  y  .  The 
corresponding  value  of  a  for  each  pulse  is  obtained  from  Equation  (3.8).  The  solid 
curves  in  Figure  5  display  pulse  speed  curves,  c  versus  a  ,  for  different  values  of 
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. . 


y  and  c  ■  0.2.  Bach  curve  has  a  well  defined  knee  (a  ,c  )  characterized  by 

v  v 

F(a  ,e,c  ,y)  ■  0  and  F  (a  ,  e,c  ,  y)  -  0  .  For  each  value  of  a  ,  rr? — r  <  a  <  a  ,  we 
v  v  c  v  v  2(1+7)  v 

have  two  pulse  solutions  with  speeds  c  (a,e,y)  and  c  (a ,  e,y)  with  c*  >  c_  .  Rinzel 

f  s  r  s 

and  Keller  (15)  showed  that  in  the  case  y  -  0  the  lower  branch,  cs  ,  corresponds  to 
unstable  pulses  while  at  the  knee  (a^c^)  the  pulse  is  neutrally  stable.  Feroe  [5) 
showed  by  numerically  evaluating  an  analytic  criterion  that  the  upper  branch,  cf  , 
corresponds  to  stable  pulses.  For  some  values  of  e  and  y  ,  other  than  those  in  Figure 
5  ,  the  pulse  speed  curve  has  no  upper  branch  or  knee;  in  such  a  case,  there  is  a  unique 
(unstable)  pulse  for  a  <  and  no  pulse  for  larger  values  of  a  . 

We  observe  that  cmin^a^'  the  minimum  value  of  c  on  each  speed  curve,  occurs  at 
a  ■  0  .  The  maximum  value  of  c  ,  cM  ,  on  each  speed  curve,  occurs  at  a  -  ^  where 

the  curve  terminates  (e.g.,  □  for  y  -  1.0)  on  the  corresponding  speed  curve  (dashed)  for 
front  solutions;  here,  z^  -  •  .  These  results  are  predicted  by  Theorem  3.1,  parts  (iii,b 
and  d). 

Fronts :  If  in  the  construction  of  a  solitary  pulse,  we  take  0  (see  Equation  (3.3)) 

then  we  will  have  a  front  as  shown  in  Figure  2A.  The  front  can  be  written  explicitly  as 


v  (z) 
c 


V 


aX^e 


- 1  4  vv  +  v , 

k-2  Pk  k 


z  <  0 


0  <  z 


where 

(3.7), 

(3.4). 


V  and  (k  =  1,2,3)  were  defined  in  (3.5),  6^  and  p^  were  defined  in 

and  the  are  the  roots  of  the  characteristic  polynomial  p(x)  defined  in 

From  Equation  (3.6)  we  find  that  the  condition  B ^  *  0  is  equavalent  to: 


p; 


(3.10)  a 
In  order  to  have  a  front  we  must  impose  the  further  restriction  that 

(3.11)  — <y  or  a  < 

1-a  f  1+7 

to  insure  that  there  are  multiple  equilibrium  states  of  system  (2.1). 
an  (3.11)  we  see  that  a  front  with  speed  c  exists  if  and  only  if  c 


By  combining  (3.10) 
satisfies 
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Pulse  (solid  curves)  and  front  (short  dashed  curves)  speed  curves 
c  versus  a,  for  different  values  of  y  and  e  =  .2.  The  long 

dashed  curve1  corresponds  to  the  rpood  curve  for  the  front  v?hen  r 


i 


'  p', 

the  corresponding  value  of  a  is  given  by  (3.10),  The  following  theorem  shows  that  a 

front  exists  for  c  sufficiently  small  or  large.  It  is  proved  in  Appendix  B  . 

Theorem  3.2?  Assume  that  the  parameters  e  and  y  are  fixed.  Then  there  exists  positive 

constants  c c2  with  c1  <  c2  such  that  if  0  <  c  <  or  c  >  c2  then  (3.12)  is 

satisfied,  and  hence  a  front  exists*  Furthermore,  lim  a  *  ^  and 

c40  2U+Y) 


if  cy  >  1 

if  cy2  ■  1 

if  ET2  <  1 


The  front  for  0  <  c  <<  1  is  unstable  (or  linearly  stable)  according  as  3a/3c|  >  0 

c~0 

(or  <  0 ) . 

We  have  solved  Equation  (3.10)  numerically  for  fixed  values  of  the  parameters 
e,  c,  and  y  •  The  short  dashed  curves  in  Figure  5  display  various  front  speed  curves, 
c  versus  a  ,  for  different  values  of  y  and  e  =  .2.  For  y  =  .25  and  y  =  1 ,  the 
speed  curves  have  two  branches.  One  expects,  as  in  the  case  of  solitary  pulses,  that  the 
faster  front  corresponds  to  a  stable  wave  while  the  slower  front  corresponds  to  an  unstable 
wave.  Notice  that  for  y  *  .25  the  speed  curve  is  not  defined  for  all  values  of  the 
speed  c  •  This  is  because  from  Equation  (3.11)  it  follows  that  a  front  cannot  exist  for 
a  >  =  .2.  These  cutoffs  for  the  speed  curve  y  =  .25  are  illustrated  by  small 

triangles. 

As  expected  from  Theorem  3.2  we  find  that  lim  a  =  -j— on  e*ch  speed  curve.  This 

is  indicated  by  the  semicircles  in  Figure  5  which  were  placed  at  the  points 

(a,c)  **  f-r-rr — ,  Ol  Note  that  „  >  0  for  y  =  .25  and  y  *  1  ,  while 

2  ( 1  +  y )  /  3c  I  c=0 

3a  I 

—l  ^  <  0  for  y  *  4  .  This  agrees  with  the  last  result  stated  in  Theorem  3.2.  The 
3c  |  c»0  ’ 

proof  of  Theorem  3.2  also  generalizes  the  results  on  limiting  behavior  for  c  small,  the 
bifurcation  of  a  slow  front  from  a  standing  front,  to  "cubic-likeM  f(v)  . 


•  v:/nr 


Figure  5  demonstrates  that  the  speed  curves  for  pulses  lie  'inside*  those  of  fronts. 
That  is,  if  the  parameters  e,y,  and  c  are  fixed  then  the  value  of  a  for  which  a  pulse 
exists  is  always  less  than  the  value  of  a  for  which  the  front  exists*  This  can  be  proven 
rigorously  by  considering  Equations  (3.8)  and  (3.10). 

The  long  dashed  curve  in  Figure  5  corresponds  to  the  speed  curve  for  the  front  when 
e  *  0  .  It  is  given  explicitly  by  c  *  ( 1-2a> [a( 1-a) ] “  ^2  .  As  the  singular  perturbation 
construction,  predicts  for  small  e  ,  the  stable  branch  of  the  speed  curves  for  the 
fronts  and  the  pulses  lie  very  close  to  the  speed  curve  for  e  »  0  . 

Backs :  It  is  not  necessary  to  present  the  analytic  construction  of  a  back  because  there  is 

a  one  to  one  correspondence  between  fronts  and  backs  under  the  change  of  variables 
(v,w)  ♦  (  ^  -  v  ,  -  w)  .  Hence  a  back  with  parameters  (a,e,c,y)  is  identified, 


under  the  above  change  of  variables,  to  a  front  with  parameters  “  a,e,c,y)  .  The 
speed  curves,  c  versus  a  with  e  and  y  held  fixed,  for  the  fronts  and  backs  are 
therefore  symmetric  with  respect  to  the  vertical  line  a  *  2 (1+y)  *  The  same  corres“ 
pondence  applies  to  pulses  and  E-pulses. 


Figure  6  offers  a  schematic  summary  of  our  numerical  and  analytic  results  concerning 
the  existence  of  known  pulse  type  and  transition  type  traveling  wave  solutions  of  the 
piecewise  linear  model.  It  shows  the  first  quadrant  of  the  y  -  c  plane  divided  into  ten 
regions.  In  each  region  there  are  four  numbers  which  indicate,  respectively,  the 
multiplicity  of  the  different  types  of  wave,  pulses,  fronts,  backs,  and  E-pulses,  which 
exist  for  those  particular  values  of  y  and  e  .  In  this  diagram  we  assume  that  the 

parameter  a  is  fixed  with  a  e  (0,  V2  )  and  we  only  consider  waves  with  positive  speed. 

a  2a 

The  constants  yQ,  y^,  and  y  are  chosen  as  before.  That  is  ,  yQ  *  — -  ,  y^  *  ^  ■  ■  , 
f  1  +  2a 

ana  *2  ■  Uii  • 

Figure  6  was  constructed  by  analyzing  various  speed  curves,  c  versus  a  and  c 
versus  y  ,  for  different  values  of  e  .  We  did  not  consider  very  large  values  of  e  . 


It  is,  of  course,  possible  that  no  waves  exist  if  e  is  too  large. 


y**Mk**'‘  * 


Whenever  just  one  front  or  back  exists  in  a  particular  reqion  it  corresponds  to  a  wave 
lying  on  the  upper  branch  of  a  speed  curve.  One,  therefore,  expects  it  to  be  stable. 
Whenever  just  one  pulse  or  E-pulse  exists  in  a  particular  region  it  corresponds  to  a  wave 
lying  on  the  lower  branch  of  a  speed  curve,  and  is,  therefore,  unstable. 

In  addition  to  the  four  types  of  waves  already  mentioned  there  also  exist  waves  with 

zero  speed,  i.e.  stationary  waves.  Note  that  these  waves  do  not  depend  on  the  parameter 

e  «  For  y  e  there  exists  a  stationary  E-pulse.  That  is  a  stationary  wave, 

V{x)  *  { v(  x)  ,u(  x)w(  x) ) ,  such  that  lim  V(x)  *  E.  For  y  >  y  there  exists  a 

I  x  | 

stationary  wave  such  that  lim  V(x)  *  0 .  This  will  be  called  a  stationary  pulse.  For 

|  x  |  -*» 

y  ■  Yj,  two  stationary  transition  waves  exist.  One  satisfies  lim  V  =  E  and 

X-f— ao 

1  *.m  V  =  0  while  the  other  satisfies  lim  V  =  0  and  lim  V  =  E.  These  represent 

X  X+-«>  X-H-w 

limiting  cases  of  the  stationary  pulse  and  stationary  E-pulse  with  infinite  pulse  width. 

Not  only  does  Figure  6  illustrate  how  many  traveling  waves  exist  for  a  given  set  of 
parameters,  but  also  how  the  waves  appear  or  disappear  through  bifurcations  as  the 
parameters  are  varied.  The  dashed  curves,  for  example,  represent  the  existence  of 
neutrally  stable  waves  out  of  which  bifurcate  the  stable  and  unstable  waves.  These 
correspond  to  the  "knees"  in  the  speed  curves  of  Figure  5.  The  dashed  curve  separating 
regions  (IX)  from  (X)  corresponds  to  neutrally  stable  fronts,  while  the  dashed  curve 
separating  regions  (I)  and  (X)  from  regions  (II)  and  (III)  corresponds  to  neutrally  stable 
pulses.  Both  curves  are  drawn  as  monotonically  increasing  funtions  of  y  because  the 
numerical  results  Indicate,  for  fixed  a  ,  that  the  minimum  value  of  y  for  which  the 
"knees"  of  the  speed  curves  in  Figure  5  exist  is  an  increasing  function  of  e  .  This  makes 
intuitive  sense  for  the  following  reason.  Suppose  the  recovery  damping  constant  y  is 
slightly  increased  from  a  critical  value  (for  which  the  pulse  is  neutrally  stable).  This 
diminishes  the  recovery  or  inhibitory  effect  of  w  so  the  system  would  exhibit  a  stable 
pulse.  Moreover,  this  pulse  would  persist  even  for  a  slightly  greater  recovery  growth  rate 
e  .  Figure  6  also  shows  that  for  fixed  e  the  neutrally  stable  pulse  occurs  at  a  larger 
value  of  y  than  the  neutrally  stable  front.  This  follows  from  a  preceding  remark  that 
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the  speed  curves  for  pulses  lie  “inside”  those  of  fronts.  It  means  that  a  stable  pulse 
occurs  simultaneously  with  a  stable  front  only  if  w  is  adequately  damped. 

The  dashed  curve  separating  regions  (IV)  and  (V)  corresponds  to  a  neutrally  stable 
E-pulse,  while  the  dashed  curve  separating  regions  (V)  and  (VI)  corresponds  to  a  neutrally 
stable  back#  Both  curves  are  drawn  as  monotonicaily  decreasing  functions  of  y  and  both 
approach  y  »  y^  as  e  ♦  0  .  This  limiting  case  has  the  speeds  of  the  E-pulse  and  back 

tending  to  zero  and  the  width  of  the  E-pulse  tending  to  infinity. 

2 

The  solid  curve  in  Figure  6  is  given  explcitly  as  ey  =  1  *  The  transition  from 
region  (VIII)  into  region  (IX)  corresponds  to  the  bifurcation  of  a  slow  unstable  E-pulse 
from  the  stationary  5-pulse.  The  transition  from  region  (VII)  into  region  (VI)  corresponds 

to  the  bifurcation  of  a  slow  unstable  pulse  from  the  stationary  pulse.  No  pulses  are  shown 
2 

for  ey  >  t  because  of  Theorem  3.1(i). 

The  other  bifurcations  occur  at  y  =  yQ  and  y  =  y^  .  As  y  increases  past  yQ 
System  (1.1)  picks  up  the  equilibrium  state  E.  For  e  sufficiently  small  we  therefore 
have  the  creation  of  the  stable  and  unstable  fronts,  a  stable  back,  and  an  unstable  E- 
pulse.  It  is,  of  course,  at  y  =  y^  that  the  discontinuity  of  f  is  reflected.  For  f 
given  by  (CUB)  the  creation  of  the  fronts,  backs,  and  E -pulses  would  not  be  so  sudden. 
Since  System  (1.1)  would  have  another  equilibrium  state  besides  0  and  E  there  may  be 
other  bounded  solutions  of  System  (1.1)  which  we  have  not  considered. 

For  y  <  y1  the  speed  of  the  stable  back  is  greater  than  that  of  the  stable  front. 
The  opposite  is  true  for  y  >  y^  .  When  y  »  y^ ,  the  stable  front  and  stable  back  have 
the  same  speed.  As  y  1  y1  the  width,  Zj ,  of  the  stable  pulse  approaches  •  In  some 

sense,  therefore,  the  stable  pulse  bifurcates  from  the  front  and  back  as  y  decreases  past 
y  =  y^.  Other  bifurcations  also  occur  at  y  =  y^  .  For  example,  as  y  t  y^  both  the 
stationary  E-pulse  and  the  slow  unstable  front  approach  a  stationary  transition  wave  which 
connects  the  rest  state  with  the  excited  state. 

Finally,  when  both  a  stable  front  and  stable  back  exist,  a  stable  pulse  also  exists 
only  if  the  back  speed  exceeds  the  front  speed;  this  necessary  condition  is  satisfied  by 


parameter  values  in  region  (III).  The  condition  is  not  sufficient  however  as  evidenced  by 
region  (X).  The  corresponding  statements  also  hold  for  E -pulses  in  regions  (IV)  and  (V) 
respectively. 
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l.  Transient  Interactions  of  Waves  and  Local  Disturbances. 

We  performed  a  number  of  numerical  experiments  to  determine  various  properties  of 
traveling  waves.  In  all  of  the  results  illustrated  here  we  used  the  piecewise  linear  model, 
f(v)  -  v  -  H( v  -  a),  with  parameter  values  a  •  .25,  e  -  .1,  and  y  -  .5;  these  para¬ 
meter  values  correspond  to  region  (III)  in  Figure  6.  For  numerical  integration,  we  used  a 
Crank-Nicholson  scheme  throughout  with  Ax  »  •  1 ,  At  *  .05,  and  a  zero  flux  condition, 

3v/3x  *  0  ,  at  both  boundaries.  For  our  choice  of  parameters  the  excited  state  is 
1  2 

and  a  stable  front,  a  faster  stable  back,  and  a  stable  pulse  exist. 

v  w  3  3 

In  describing  the  numerical  experiments  it  will  be  convenient  to  refer  to  different 
phases  of  the  traveling  waves.  We  shall  refer  to  the  sharp  transition  from  the  left  branch 
to  the  right  branch  of  f(v)  for  the  front  or  pulse  as  "the  leading  edge",  and  the  sharp 
transition  from  the  right  branch  to  the  left  branch  of  f(v)  for  the  back  or  pulse  as  the 
"trailing  edge".  These,  of  course,  correspond  to  the  inner  regions  in  the  singular 
perturbation  construction. 

Formation  of  a  Pulse  from  Square  Initial  Data.  The  formation  of  a  pulse  from  specified 

initial  conditions  is  illustrated  in  Figure  7  .  Spatial  profiles  v(x,t)  versus  x  (Figure 

7A)  and  w(x,t)  versus  x  (Figure  7b)  are  shown  for  discrete  t:  t  ■  5p  #  0  <  p  <  8  . 

For  initial  data  we  used  square  steps  of  unit  height:  v(x,0)  =  Xrc  ici  and  w(x,0) 

15 , 35  J 

=  Xm  c  t  ci  where  \f  ..  is  the  characteristic  step  function  of  the  interval  fa,b]. 
[2.5,7*5J  La,bJ 

Note  the  asymmetrical  distribution  of  v(x,0),  w(x,0)  with  respect  to  each  other;  the 
step  in  w(x,0)  overhangs  that  in  v(x,0).  This  helps  ensure  the  formation  of  a  pulse 
rather  than  two  fronts  traveling  in  opposite  directions. 

We  see  that  the  leading  edge  in  the  front  part  of  the  pulse  is  clearly  formed  at 
t  *  t^  and  proceeds  to  travel  with  constant  velocity.  The  trailing  edge  in  the  back  part 
of  the  pulse  also  forms  quickly  but  begins  by  traveling  very  slowly  (tQ  <  t  <  t2).  The 
trailing  edge  then  accelerates  until  it  is  actually  moving  faster  than  the  leading  edge 
(t1  <  t  <  t^).  The  trailing  edge  then  slows  down  and  its  speed  approaches  that  of  the 
leading  edge.  Then  the  leading  and  trailing  edges  travel  with  the  same  speed 
( t^  <  t  <  tg)  and  a  pulse  is  formed. 
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Formation  of  a  pulse  from  square  initial  data.  Spatial  profiles 
v(x,t)  versus  x  are  shown  for  discrete  t  :  t  -  5p,  0  <  p  <  I 


Figure  7B.  Formation  of  a  pulse  from  square  initial  data.  Spatial  profiles 
w(x,t)  versus  x  are  shown  for  discrete  t  s  t  *  5p,  0  <  p  <  8. 


In  order  to  understand  why  the  pulse  is  formed  in  this  way  recall  that  in  the  singular 


perturbation  construction  of  a  traveling  wave  the  recovery  variable,  w  ,  is  assumed  to  be 
constant  in  the  inner  region.  The  speed  of  the  leading  edge  is  a  decreasing  function  of 
this  constant  while  the  speed  of  the  trailing  edge  is  an  increasing  function  of  this 
constant.  Ahead  of  the  pulse  w  is  always  very  close  to  zero,  and  hence  the  leading  edge 
propagates  with  constant  speed.  The  trailing  edge,  however,  does  not  see  a  constant  value 
of  w  ,  and  therefore  its  speed  varies. 

In  the  region  between  the  trailing  and  leading  edges  the  spatial  derivatives  of  v 
are  very  small,  and  hence  the  solution  will  follow  very  closely  space  clamped  dynamics. 

That  is,  the  solution  (v,w)  will  evolve  according  to  the  ordinary  differential  equation: 

v  *  f(v)  -  w 
w t  ■  e(v  -  yw) 

with  initial  conditions  v(0)  =  1  and  w(0)  »  0  .  It  follows  that  in  the  middle  region 
the  solution  <v,w)  becomes  very  close  to  the  excited  state  (Ey,  Ew) . 

The  velocity  of  the  trailing  edge  initially  increases  until  this  edge  sees  a  value 

of  w  very  close  to  w  *  ?w.  Because  of  our  choice  of  parameters  the  trailing  edge  is  now 

traveling  faster  than  the  leading  edge.  After  the  trailing  edge  passes  x  3  35  its  speed 
decreases  because  then  it  advances  into  a  region  where  w  is  decreasing.  The  pulse  is 
formed  when  the  trailing  edge  slows  to  the  same  speed  as  the  leading  edge.  Then,  be-cause 
both  edges  have  the  same  speed,  the  trailing  one  always  sees  a  constant  value  of  w. 

For  a  singular  perturbation  construction  of  how  traveling  waves  are  formed  from 
initial  data  see  Keener  [9]  . 

Collision  of  a  Pulse  and  a  Front.  In  Figure  8  we  take  for  initial  data  (dashed  curves)  a 
pulse  moving  to  the  right  and  a  front  moving  to  the  left.  Figure  8A  shows  profiles,  v(x,t) 
versus  x  ,  and  Figure  8B  shows  profiles  w(x,t)  versus  x  ,  for  discrete  t;  t  *  5p, 

0  <  p  <  7  .  When  the  two  leading  edges  reach  each  other  the  dynamics  are  very  similar  to 

the  collision  of  two  fronts.  This  brings  the  region  near  the  fronts  to  the  excited 


-29- 


X 


state  (t^  <  t  <  t^).  When  the  trailing  edge  of  the  pulse  enters  this  excited  region  it 

increases  in  velocity  and  evolves  into  a  back  moving  to  the  right  (t^  <  t  <  t?).  The 

reason  that  the  trailing  edge  increases  in  velocity  is  that  the  value  of  w  at  the  excited 

state  is  greater  than  the  value  of  w  at  the  trailing  edge  of  a  pulse.  The  singular 

perturbation  construction  showed  that  the  velocity  of  the  trailing  edge  of  the  pulse  is  an 

increasing  function  of  the  value  of  w  at  which  the  trailing  edge  takes  place. 

Reflections  of  a  Traveling  Wave.  Suppose  that  (a(x,t),  0(x,t))  is  the  solution  of  System 

(1.1)  with  initial  data  (r^(x),  ^(x))*  We  define  (n^(x),  r^(x))  to  be  a  local 

disturbance  if  lim  {sup  (|a(x,t)|  +  l8(x,t)))  *  0.  This  is  a  local  disturbance  with 
t+«  xen 

respect  to  the  steady  state  (v,w)  »  (0,0);  one  may  also  define  analogously  a  local 
disturbance  with  respect  to  (Ev,Ew).  We  say  (v(x,t),  w(x,t))  is  a  reflection  of  the 
traveling  wave  solution  (v(x,t),  w(x,t)),  a  reflection  caused  by  a  local  disturbance,  if 

(i)  For  some  time  T  >  0  and  some  local  disturbance  (n^(x),  n2(x))» 

(v(x,t),w(x,t) )  is  the  solution  of  systems  (1.1)  with  initial  data 
( v(  x,T)  +  n^x),  w(x,T)  +  n  ( x) ) ,  and 

(ii)  There  exists  a  traveling  wave  solution  (Vj(x,t),  w^(x,t)  moving  in  the 

opposite  direction  of  (v,w)  such  that 

lim  (sup(|v(x,t)  -  Vj(x,t)|  +  I w( x , t )  -  wt(x,t)|)}  *  0  . 
t  xGR 

We  gay  that  it  is  possible  to  reflect  the  traveling  wave  (v,w)  if  there  exists  a 
reflection  of  (v,w). 

Our  numerical  experiments  indicate  that  for  given  values  of  the  parameters 
a,  e,  and  y  it  may  be  possible  to  reflect  some  traveling  waves  but  not  others.  Figure  9 
shows  a  front  being  reflected  into  a  back  with  profiles  v(x,t)  versus  x  (Figure  9A)  and 
profiles  w( x , t )  versus  x  (Figure  9B)  for  discrete  t:  tp  *  2.5p  ,  0  <  p  <  17,  The 

first  three  profiles  show  a  front  moving  to  the  right.  At  t  *  t3  (dashed  curves)  we 
increased  the  values  of  w  over  a  short  interval  just  ahead  of  the  leading  edge.  This 
certainly  represents  a  local  disturbance.  In  fact  one  sees  that  in  the  disturbed  region 
the  solution  slowly  returns  to  rest.  As  the  leading  edge  enters  the  region  of  increased 
w  it  glows  down  (t^  <  t  <  t10)  and  then  collapses  to  become  a  back  moving  towards  the 


left  (t1Q  <  t  <  t17).  Notice  that  the  back  travels  with  a  speed  greater  than  the  front,  as 
indeed  it  should  for  parameter  values  from  region  (III)  of  Figure  6.  This  is  also 
consistent  with  the  singular  perturbation  prediction  which  we  note  gives  reasonably 


accurate  estimates  for  the  back  and  front  speed  if  e  is  small. 

It  may  be  surprising  to  discover  that  under  the  corresponding  circumstances  a  back 
will  not  reflect  to  become  a  front.  This  is  shown  in  Figure  10  with  Figure  10A  showing 
profiles  v(x,t)  versus  x  and  Figure  10B  showing  profiles  w(x,t)  versus  x  for 
discrete  t:  t  ■  2.5p,  0  <  p  <  14.  The  first  two  profiles  show  a  back  moving  towards  the 
right.  At  t  -  tj  (dashed  lines)  we  sharply  decreased  w  over  an  interval  directly 
ahead  of  the  trailing  edge  of  the  back.  This  is  certainly  a  local  disturbance  (with 
respect  to  (Ev,  Ew)).  In  one  sees  that  in  the  disturbed  region  the  solution  slowly 

returns  to  the  excited  state.  Upon  entering  the  disturbed  region  the  trailing  edge  slows 
down  and  then  apparently  stops,  waiting  for  the  disturbance  ahead  of  it  to  return  to  the 
excited  state  ( t^  <  t  <  t^).  The  back  then  continues  to  move  in  its  original  direction 
(t9  <  t  <  t15). 

Our  numerical  experiments  suggest  that  if  the  parameters  a,  e,  and  y  lie  in  a 
certain  subset  of  values  for  which  the  back  speed  exceeds  the  front  speed  then  it  is 
possible  to  reflect  a  front  but  impossible  to  reflect  a  back  by  local  disturbances.  To 
understand  this  behavior  we  formulate  and  study  an  idealized  model  problem.  We  consider 
first  the  case  of  a  back  moving  rightward  toward  a  local  disturbance  of  decreased  w 
(Figure  10).  As  one  recalls  from  the  singular  perturbation  construction,  the  back  speed  is 
determined  by  the  value  of  w  directly  ahead  of  it  and,  moreover,  this  speed  decreases 
with  decreasing  w  .  Hence  one  expects  intuitively  that  if  w  falls  abruptly  in  the 
disturbed  region  then  the  back  will  decelerate  suddenly  and  virtually  halt.  An  appropriate 
formal  perturbation  calculation  for  this  transient  stage,  which  includes  the  deceleration, 
temporary  halt,  and  initiation  of  either  reflected  or  continued  propagation,  is  yet  to  be 
worked  out.  Our  approach  is  to  observe  from  the  numerical  results  during  this  transient 
state,  T  <  t  <  Tj  say,  that  v  and  w  at  the  "halting  location”  Xq  have  values  close 
to  Ev  and  Ew  .  Thus  we  consider  the  model  problem  in  which  v,w  satisfy  the  equations 
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(4.  la) 


V  »  V  +  f  ( v)  -  w 
t  XX 

wfc  -  e(v  -  yw)  for  x  <  xQ  ,  t  e  (T ,  T^  ) 

with  the  boundary  conditions; 

(4.1b)  (vfx^t),  w(xQft))  -  <Ey,  Ew) 

(4.1c)  (v,w)  -►  (0,0)  as  x  ♦  -• 

and  with  the  approaching  back  profile  as  initial  conditions. 

We  remark  that  one  can  also  motivate  the  boundary  condition  (4.1b)  as  follows. 

Let  V  denote  the  disturbed  region  in  Figure  10.  Shortly  after  the  disturbance  occurs, 
v  >  Ev  ,  w<Ew  in  P  and  so  v  diffuses  outward  across  3P  .  Hence,  just  outside  P  , 
wfc  >  0  and  one  observes  a  local  increase  in  w  there.  The  local  spread  of  v  also  means 
that  wfc  is  less  positive  just  inside  3P  than  it  is  toward  the  middle  of  P.  This 
explains  the  under/overshoots  in  the  w-profile  (Figure  10B)  that  one  sees  in  crossing  3P 
outward.  From  this  it  also  follows  that  between  the  under  and  overshoot,  w  *  E 

w 

(at  x  *  xQ,  as  in  the  preceding  paragraph) ,  at  least  for  a  while.  Finally,  if  the  back 
approaches  soon  then  v  >  E^  just  inside  and  v  <  just  to  the  left  of  3P  so  that 

v  *  E^  near  xQ  .  The  above  argument  does  not  apply  if  the  local  disturbance  is  not 
fairly  sharp.  But  of  course,  one  can  probably  not  induce  reflections  by  disturbances  which 
are  not  sharp.  At  the  same  time,  the  singular  perturbation  treatment  would  not  predict 
halting  if  the  local  disturbance  were  slowly  varying  in  x  ;  ir  such  a  case,  the  wave  could 
decelerate  sufficiently  without  stopping  as  the  disturbance  wears  off. 

One  expects  that  the  temporary  halt  of  a  back  may  be  viewed  as  the  approach  of  the 
solution  of  (4.1)  to  a  time-independent  steady  state  solution.  The  steady  state  solution 
(q(x),  p(  x )  )  must  satisfy 

q"  +  f(q)  -  p  *  0 

q  -  yp  -  0 

with  boundary  conditions 

(q(  Xq  )  ,p(  Xg  ) )  -  (Ey,  Ew)  and  (q( -®)  ,p( -«) )  *  (0,0). 
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q(x)  which  satisfies  the  boundary 


Since  p(x)  »  q(x)/y  in  {-*»,  xQ)  we  seek  a  function 
value  problem: 


(4.2) 


q"  +  f<q)  -  q/y  -  0  in  (-®,x0) 


q<-°°)  =  0  ;  q(xQ)  =  . 

For  y  e  (y^y^  problem  (4.2)  has  a  unique  solution.  This  can  be  seen  by  considering  the 
phase  plane  portrait  for  the  equation  q"  +  f(q)  -  q/y  =0  as  shown  in  Figure  11. 

Finally,  the  steady  state  will  be  approached  only  if  it  is  stable  as  a  solution  to  the 
System  (4.1).  In  Theorem  4.1  we  demonstrate  (see  Appendix  C)  the  stability  or  instability 
of  (q,p). 

Theorem  4.1. 

2 

Suppose  yQ  <  ^  <  '  or  equivalently  1  <  2a /Ev  <  2  .  If  ey  >  1  ,  then  (q(x), 

2 

p(x))  is  linearly  stable  as  a  solution  to  (4.1).  For  ey  <  1  ,  (q(x),  p(x))  is 

linearly  stable  when 


(4.3)  2a/Ev  >  1.278... 

but  unstable,  for  given  e,y,  if 


(4.4)  0  <  2a/Ev  -  1  <<  1  . 

We  conjecture  that  when  the  back  is  faster  than  the  front  and  (q(x),  p(x))  is  stable 
then  the  back  cannot  be  reflected  by  local  disturbances. 

A  few  comments  on  Theorem  4.1  are  in  order.  The  sufficient  condition  (4.4)  which 
guarantees  instabilitv  is  an  asymptotic  condition.  Precise  conditions  are  those  for  which 
inequality  (C.11)  of  Appendix  C  is  satisfied.  As  the  proof  demonstrates,  the  loss  or  gain 
of  stability  as  parameters  are  varied  leads  to  Hopf  bifurcation  of  periodic  solutions. 

This  applies  to  more  general  and  related  problems  for  bistable  systems  which  we  will  pursue 

2 

in  further  work.  Finally,  note  that  the  statement  in  Theorem  4.1  for  ey  >  1  holds  for 
any  "cubic- like"  f(v)  ;  see  Appendix  C  • 

Our  computed  examples  are  for  parameter  values  which  satisfy  inequality  (4.3).  Hence, 
corresponding  to  our  model  problem  (4.1),  we  identify  the  halted  back  in  Figure  10  with  the 
steady  state  solution  (q,p)  of  (4.1).  The  back  waits  at  this  steady  state  for  the  local 


.gure 


disturbance  to  decay  sufficiently.  When  this  has  happened  the  back  can  proceed  through  the 


disturbed  region.  For  parameter  values  such  that  the  steady  state  (q(x),p(x))  exists  but 
is  unstable  we  would  not  predict  that  it  is  impossible  to  reflect  the  back. 

For  the  analogous  problem  (Figure  9),  the  front  does  reflect,  rather  than  stop  and 
wait  for  the  local  disturbance  to  decay,  because  the  corresponding  steady  state  solution 
does  not  exist.  To  be  more  precise,  the  model  problem,  in  this  case,  is  the  system  of 
Equations  (4.1a)  with  the  boundary  conditions 
(4.5a)  (v(xQ,t),  w(xQ,t))  -  (0,0)  , 

(4.5b)  ( v, w)  ♦  (E  ,E  )  as  x  ♦ 

V  w 

rather  than  (4.1b,c).  We  would  seek  a  steady  state  q(x),  with  p(x)  -  q(x)/y  ,  as  a 
solution  to  the  boundary  value  problem 


(4.6) 


q"+  f ( q)  -  q/Y  -  0 


q(-<»)  =  Ev  ,  q<XQ>  3  0  • 


From  the  phase  plane  portrait  ( Figure  11),  we  see  that  this  problem  does  not  have  a 
solution. 

I<et  us  examine  further  the  transient  phase  of  the  reflection  in  Figure  9.  Observe 
that  the  leading  edge  is  stopped  by  the  local  disturbance  at  t  *  t^  and  does  not  move 
for  t4  <  t  <  t1Q.  It  cannot  move  to  the  right  because  of  the  local  disturbance  in  w  . 

It  also  cannot  move  to  the  left.  If  it  did  it  would  resemble  the  trailing  edge  of  a  back 
moving  to  the  left.  Since  the  value  of  w  immediately  to  the  left  of  this  trailing  edge 
is  very  close  to  zero  this  would  be  impossible.  The  leading  edge  of  the  front  must 
therefore  wait  at  x  **  x^  for  either  the  values  of  w  directly  ahead  of  it  to  decrease 
sufficiently  so  that  it  can  move  to  the  right,  or  for  the  values  of  w  directly  behind  it 
to  increase  sufficiently  so  that  it  can  move  to  the  left.  In  Figure  9b  we  see  that  the 
value  of  w  steadily  increases  behind  the  leading  edge  for  <  t  <  tg.  When  these 

values  of  w  are  large  enough  the  leading  edge  collapses  to  form  a  back  moving  to  the 
left.  Hence,  for  a  front  to  reflect  the  local  disturbance  must  be  sufficiently  large  so 
that  the  leading  edge  is  stopped  long  enough  for  the  values  of  w  behind  it  to  increase 
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sufficiently.  Note  that  the  values  of  w  behind  the  leading  edge  would  be  unable  to 
increase  sufficiently  if  there  existed  a  steady  solution  of  (4.5). 

The  analysis  we  have  presented  so  far  holds  for  ye  (Yq/Y-j)*  Only  for  these  values 
of  y  does  the  back  speed  exceed  the  front  speed.  On  the  other  hand#  the  front  is  faster 
than  the  back  if  ye  (y^y^)*  In  this  case  there  exists  a  solution  of  the  boundary  value 
problem  (4.6)  but  there  does  not  exist  a  solution  of  Equation  (4.2).  Hence  we  conjecture 
that,  for  a  subset  of  parameter  values  for  which  the  front  speed  exceeds  the  back  speed,  it 
is  always  possible  to  reflect  a  back  but  not  a  front. 

These  results  indicate  that  for  y  <  y^ ,  the  rest  state  in  some  sense  dominates  the 
excited  state.  For  y  <  y^  the  only  type  of  stable  traveling  wave  which  brings  the  system 
to  the  excited  state  is  a  front  which  can  be  reflected  by  means  of  a  local  disturbance.  For 
a  subset  of  this  range  of  parameters,  when  (q(x),  p(x))  is  stable,  the  back  and  pulse 
cannot  be  reflected.  Hence,  if  the  system  is  being  switched,  as  t  ♦  •  ,  towards  the 
excited  state  then  it  is  possible,  by  means  of  a  local  disturbance,  to  return  the  system  to 
rest.  However,  if  the  system  is  switching  towards  the  rest  state  then  it  is  impossible,  by 
means  of  a  local  disturbance,  to  bring  the  system  to  the  excited  state.  For  y  >  y  the 
corresponding  opposite  statement  is  true  and  the  excited  state  is  dominant. 

As  mentioned  in  the  introduction,  our  consideration  of  reflections  by  local  dis¬ 
turbances  also  leads  to  expectations  about  the  outcome  of  a  front  or  back  which  approaches 
a  boundary  where  v,w  are  held  fixed  at  either  (0,0)  or  (Ey,  Ew)  •  This  is  because  our 
model  problems  (4.  1a,b,c)  and  (4.1a,  4.5a, b)  incorporate  Dirichlet  boundary  conditions 
explicitly.  Suppose  parameters  are  such  that  the  back  speed  exceeds  the  front  speed. 

Then,  because  (4.6)  has  no  solution,  we  would  expect  a  front  to  be  reflected  from  a 

boundary  where  (v,w)  *  (0,0).  Next  consider  a  back  approaching  a  boundary  where 

(v,w)  *  (Ey,  6^).  If  the  steady  state  solution  (q(x),  p(x))  to  (4.1)  is  stable,  then  we 

expect  the  back  will  tend  to  (q(x),  p( x) )  and  not  be  reflected.  On  the  other  hand,  if 

2 

(q(x),  p ( x ) )  is  unstable  ( ey  <1  and  inequality  (4.11),  cf.  Theorem  4.1)  then  we  expect 
that  the  back  might  be  reflected.  Indeed,  we  have  obtained  numerical  solutions  of  the 
model  problems,  for  parameter  values  in  this  latter  regime,  which  illustrate  that 
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reflection  may  take  place  for  a  back  approaching  an  excited  boundary  as  well  as  for  a  front 
approaching  a  resting  boundary. 
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Appendix  A  -  Proof  of  Theorem  3.1. 

Throughout  this  discussion  we  assume  that  the  parameters  e  and  y  are  fixed  and  set 
h(s,c)  *  F(a,e,c,y).  We  determine  for  which  values  of  c  there  exists  a  root,  s  ,  of  the 
equation 

h(s,c)  -  0 


such  that  s  e  (0,1). 
Note  that. 


(A.1)  h( 1 ,c)  *  —  <1,c)  *  0 

as 

for  all  values  of  c  .  Furthermore,  for  any  fixed  value  of  c  there  can  exist  at  most  one 
3^h  *  * 

root  of  the  equation  — r  (s  ,c)  *  0  for  a  >  0  .  If  a  root  exists  then  it  is  given  by 
3s2 


* 


s 


p'3 


V(VV 


In  the  case  we  must  choose  a  branch  of  the  logarithm  to  determine  s  .  This, 

along  yith  (A.1),  implies  that  for  a  given  value  of  c  at  most  one  pulse  can  exist  for  a 
uniquely  determined  value  of  s  and,  therefore,  from  (3.8),  for  a  uniquely  determined 
value  of  a  .  Furthermore, 

Lemma  A. 1 .  Assume  that  the  parameters  e, y  ,  and  c  are  fixed.  Then  a  necessary  and 
sufficient  condition  for  pulses  to  exist  is  that  either 

,2 

(a)  h(  0 ,  c)  >  0  and  ^h(1,c)  <  0 

3s 

(A. 2)  or 

a2 

(b)  h( 0, c)  <  0  and  h(1,c)  >  0. 

3s2 

In  Lemmas  A. 3  and  A. 4  we  determine  for  which  values  of  c  the  functions  h(0,c)  and 

£ 

— -  h(1,c)  are  either  positive  or  negative.  We  first  prove  some  properties  of  a  ,  the 
3s  1 

positive  root  of  the  characteristic  polynomial  p(x)  defined  in  (3.4).  Recall  that 

2  1  +  y 

O  i  -  . 

Y 

Lemma  A. 2.  Assume  that  the  parameters  e  and  y  are  fixed.  Then 

(a)  lim  a.  =  a 

c+0 
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(b) 


s_l 

3c  1  c*0 

2e,2 

(c) 

>  c  for  all 

c  >  0 

(d) 

2 

if  ty  >  1  then  >  o 

for  all 

(e) 


if 


ey2  <  1 


then,  ol 


<  o  for  0  <  c  <  a 


lZ£L 


/y (1+CT) 


*  o  for  c  ■  a 
>  o  for  c  >  a  . 

Proof.  Parts  (a)  and  (b)  follow  from  substituting  an  expansion  of  in  powers  of  c 

into  the  characteristic  polynomial  (3.4).  For  (c)  note  that  p(c)  *  -c  -  <  q  . 

c 

Since  p(x)  >  0  for  x  >  the  result  follows. 

We  prove  (d)  by  showing  that  p(c)  <  0  for  all  c  >  0  .  Now,  p(  c)  * 

Vo  2  A  1-£Y^  2 

a  2  (o  -  co  -  ( 1+ey) )  •  Hence,  p(o)  -  0  if  c  =■  o  - - 1  -  —  .  If  ey  >  1  then 

/y( 1+cy) 

p(o)  is  never  equal  to  zero  for  any  positive  value  of  c  .  Therefore,  either  p(  o)  <0 

or  p(  a)  >0  for  all  c  >  0  .  It  follows  that  either  >  a  or  <  a  for  all  c  >  0. 

But,  by  part  (e)  >  c  for  all  c  •  Choosing  c  >  a  the  result  follows. 

2 

To  prove  (e)  assume  that  ey  <  1 •  Hence  o  >  0  •  By  inspection, 
p(a)  <  0  for  c  >  a  .  Hence,  a  <  at^  because  p(x)  >0  for  any  x  >  .  Similarly, 

p(o)  >  0  for  c  <  o  so,  in  this  case,  o  >  •  // 

We  now  study  the  function  h(0,c)  which  can  be  written  explicitly  as: 


(A. 3) 


h(0,c) 


2  - 


First  note  that  the  condition  h<0,c)  >  0  is  equivalent  to 

2  2  PY^  i 

(A. 4)  H(c)  -  (c  -  e-r)a*  +  (  —  )  2ca +  (e  +  ey)  <  0  . 

2 

(A. 4)  is  obtained  by  substituting  into  (A. 3)  the  explicit  formulas  for  p^  ,  o  ,  and  6^ 
and  then  using  the  fact  that  p(  )  *  0  • 

Lemma  A. 3;  Assume  that  the  parameters  t  and  y  are  fixed. 

-  2 

(a)  There  exists  a  positive  constant  c^  such  that  if  ty  <1  and  c  e  (0,^)  then 
h(0,c)  >  0  . 
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(b)  There  exists  a  positive  constant  c such  that  if  ey  <  1  and  c  >  then 
h(0,c)  <  0  . 

(c)  If  ey2  >  1  then  h(0,c)  <  0  for  all  c  >  0  . 


From  Lemma  A. 2 (a)  and  (b)  it  follows  that  H(0)  -  0  and  H*(0)  -  f  ^  —  )<  0. 

1 v  y  * 

Part  (a)  now  follows  immediately,  (b)  is  true  because  ^  c  for  c  >>  1  so  that 
lim  H{c)  ■  +•  .  We  prove  (c)  by  showing  that  H(c)  >  0  for  all  c  >  0  •  From 
p(o^)  •  0  ,  we  replace  the  first  and  third  terms  of  H(c)  by  c  aj  •  c(ey+1)a^.  Now 


H(c)  may  be  written  as 


H(c)  -  co1  (a2  -  o2  +  ) 


Hence,  from  Lemma  A. 2(d),  H(c)  >  0  .  // 

32  a2  pi 

We  now  study  — -  h(1,c)  which  can  be  written  explicitly  as  — r  h(1,c)  -  — —  -  2. 

3s  32  3s2  a261 

Lemma  A. 4.  A  necessary  and  sufficient  condition  for  — —  h(1,c)  <0  is  that 

3s 


v2  +  /XZJ7T 


(cy+e)  1  <  c(  ey+c)  *  +  /( c  +  ey)  +  c  +  ey  • 


Proof :  If  we  substitute  the  explicit  expressions  for  p^  and  61  into  that  of  — r  h(1»c) 

2  3sZ 

and  use  the  fact  that  p(  a,  )  *  0  we  find  that  — -  h(1,c)  <0  is  equivalent  to 


which  is  equivalent  to 


This  can  be  written  as 


i£U^2  V2< 

ey  +  c 

ey  +  c 


r  -  2  c  r  -  ut+u  <  o 


r-  )V 


Finally,  (A.7)  is  equivalent  to 


ey  +•  c 


(  )1/2  -  r  <  c  ♦  /c2  ♦  «Y  ♦  i 


ey  +  c 


( ey  +  e )  ^2  <  c(  ey  +  c2 )  ^2  +  /(c2  +  ey)2  +  c2  +  ey  . 


3  3 

h(1,c  ,  )  ■  -r*  h(1,c  .  )  *  — r  h(1,c  .  )  *  0  •  Therefore,  Ms,c  .  i  has  a  triple  root 
win  3s  min  ~  2  min  m*n 

oS 

at  s*1  and  the  result  follows. 

We  now  wish  to  define  cmax  .  Recall  that  h<0,cmin)  >  0  and,  from  A. 3(b), 

h(0,c)  <  0  for  c  >  .  Let  c  -  inf {c  >  c  .  ,  h(0,c)  *  0}.  Our  previous  discussion 

2  max  min 

shows  that  pulses  exist  for  c  e  ( cm^n#  cmax^ “  Note  that  the  condition  h(°'cmax)  *  0 

pi  2 

implies  that  —  -  2o  at  .  Furthermore,  from  (3.8),  it  follows  that  at  c  -  cmax/ 


0*s 


a 


2q2  2(1*t> 


That  is,  lim  a  *  'zrrr* — r  •  Since  s  *  e  we  also  conclude  that  lim  z„  -  •  . 

etc  2(1+Y)  etc  1 

max  max 


It  remains  to  prove  that  as  ey  t  1  both  c  in  and  cM  I  0  . 


Recall  that 


11(0,0. _ )  -  0  •  Therefore,  H(c,.  )  *  0  where  H(c)  has  been  defined  in  formula  (A.4). 

max  max 

2  2 

Moreover,  at  ey  »  1  ,  H(0)*HM0)*0  and,  for  e*y  >  1#  the  proof  of  Lemma  A.  3(c) 


shows  that  H(c)  >  0  for  all  c  >  0  «  Hence,  cmax  "  0  must  be  a  double  root  of  the 

2 

equation  H(c)  =*  0  at  *1  and  the  result  follows. 
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Assume  that  the  parameters  e  and  y  are  fixed.  Recall  that  a  necessary  and 
sufficient  condition  for  a  front  to  exist  is  that  the  speed,  c  ,  satisfies  the  inequality 


12 
7~  >  o 


p; 


From  the  proof  of  Lemma  A. 3(a),  we  conclude  that  lim  h(0,c>  -  lim  a  - 


0  .  That 


c+0 


c+0  a  6 


2  2 

is,  lim  — —  -  2a  >  a  .  Therefore,  a  front  exists  for  c  sufficiently  small. 
c+0  61 


P1  2  2 

Furthermore,  if  c  is  sufficiently  large,  then  -j-  >  2a  >  a  and  a  front  exists.  This 

is  an  immediate  consequence  of  Lemma  A. 3(b)  and  (c). 

We  next  show  that  lim  a  =  — *  -J,  .  in  the  preceding  paragraph  we  saw  that 

c+0  2c  211  V 


P\ 

lim  -r± 
c+0  61 


2c  .  The  result  now  follows  from  Equation  (3.10). 
It  remains  to  prove  that 


>  0 

if 

2 

ey 

<  1 

(B.1) 

3aj 

8c|c”0 

-  0 

if 

2 

ey 

-  1 

<  0 

if 

2 

cy 

>  1  . 

It  will 

be  easier  to 

consider 

c  as  a 

function 

of 

y  • 

Note,  for  PWL,  that  at  < 

Y  » 

*  which 

1-2a 

implies 

that  3y/9a  |  > 

C”U 

0  . 

Hence,  (B.1)  is  equivalent  to 

>  0 

if 

2 

ETfl 

<  1 

(B.2) 

ill 

3c| c*0 

*  0 

if 

2 

EY1 

*  1 

<  0 

if 

2 

EY1 

>  1  . 

The  result  (B.2)  can  be  verified  algebraically  for  PWL.  In  the  following  proof  we 
require  only  that  f  is  "cubic-like"  as  shown  in  Figure  4.  For  PWL  and  CUB,  y^  is 
defined  in  Section  2.  More  generally,  when  f  may  not  have  the  symmetry  properties  of  PWL 


or  CUB,  then  y  is  the  value  of  y  such  that 


Ev 

0  -  /  (f(v)  -  v/y^dv. 


In  this  more  general  case,  one  should  remember  that  Sa/dcl^g  an<*  Oy/Se I wm  have 

opposite  sign  if  3y/3a|c_Q  <  0  * 

We  first  write  y,  v  ,  w  and  f(v  )  as  power  series  in  c 
c  c  c 

Y<C)  -  Y,  ♦  CY,  +  c2y2  +  ••• 

V  (z)  «  ^n(z)  +  C*>  (z)  +  ••• 

C  vl  1 

( B.3 ) 


W  (z)  -  +  c4>,  <*>  +  ••• 

C  1/  1 

f(vc>  *  f(v?Q)  +  +  •••  • 


Recall  that  (vc*wc)  solves 


(B.4) 


crv •  -  v"  +  f(v  )  -  w 
c  c  c  c 


cw'  »  e(v  -  yw  ) 
c  c  c 


with  boundary  conditions 

(v  (-«),  w  (-»))  *  (0,0) 
c  c 

(v  (+*),  w  (+•))  -  ( E  , E  )  . 
c  C  V  w 

Substituting  (B.3)  into  (B.4)  and  equating  like  order  terms  we  find  to  the  zeroth  order 
that: 

0  -  +  f<*0)  -  *0 

.  0  -  £<*0  - 

(V>0<—  >/♦„<—>>  -  (0.0)  and  (V>0(+»),i|.0(+«'))  -  (EV’EW>‘ 

Hence,  4»0<z)  *  *o<X)/Y1  .  Therefore,  *Q(z)  is  a  solution  of  the  ordinary  differential 
equation 

0  «  +  f(*Q)  -  V>Q 
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1 


with  boundary  conditions 


*  <— >  *  o  ,  *>0(+»)  *  Ev  . 


It  is  not  hard  to  show  that  such  a  solution  exists  and  is  a  monotonically  increasing 
function. 

We  now  equate  rirst  order  terms  to  obtain 


( B.5 ) 


-  f,<v*i  -  *1 

*0  ■  E<*1  '  (Vl  +  Vo** 


This  implies  that 


-  Vo  *0 

1  “  ’  E^i 


and,  hence,  ^  is  a  solution  of  the  equation 


( B.6 ) 


"1 


1  v  '1 


+  f,(V*i  -  -  *5  (1  -  —  >  -  -f  *o 

1  EY1 


Since  is  a  solution  of  the  homogeneous  equation: 


\p ' 

<^)"  +  <^)’  f'(^0)  -  Y  -  0  , 


a  necessary  and  sufficient  condition  for  (B.6)  to  have  a  solution,  fy  ,  is  that 


( B.  7 ) 


•  ^o(i  -  -V  >  -  -?  v  - 0  • 


Here  <u,v>  =  f*  u  v  dx.  How  ,  vpQ)  =  ~  .  Therefore,  (B.7)  becomes 


(1  -7 =  0 

£Yi  2Y, 


! 

i 

! 


or 
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2[l  - 


eY, 


'  V 


from  which  the  result  follows. 

We  note  that  the  same  proof  applies  for  the  emergence  of  the  slow  back  solution  from 

_  2 
c  *  0  except,  in  that  case,  y^  takes  the  opposite  sign  because  «  -E^/2  • 

To  treat  stability  we  introduce  a  traveling  coordinate  frame  (z,t),  z  »  x  +  ct  ,  and 

seek  exponential  solutions  (e^  V(z),  e^Wfz))  to  the  linear  variational  equation. 

Hence,  X,  V,  W,  satisfy 

XV  =  V"  -  cV’  +  f * ( v  ( z) ) V  -  W 
c 

(B.8) 

XW  =  -cW  +  e(V-  yw)  . 


As  in  the  proof  above,  we  follow  a  perturbation  argument  for  small  c  and  write 

V(z)  *  VQ(z)  +  c  V^z)  +  ••• 


W(z)  =  WQ(z)  +  c  W^z)  + 

X  “  XQ  +  c  Xf  +  ••• 


To  lowest  order  we  find  for  y  »  y^  +  0(c)  that  W^ 


xo+E^  0  ' 


and 


(B.9)  VJj  +  (u0  +  f*(^0<z))l  VQ  =  0  ,  z  G  R+ 

where 

(B.10)  W0  -  -Xg  *  e/(  Xg  +  er,)  . 

The  perturbation  procedure  is  straightforward:  X^  is  obtained  from  the  solvability 
condition  for  the  next  order  equation,  etc.  Therefore  we  concentrate  on  (B.9),  a 
Schroedinger  eigenvalue  problem.  If  (B.9)  has  a  bounded  solution  with  Re  XQ  >  0  then 
(vc,wc)  is  unstable  for  0  <  c  <<  1.  On  the  other  hand,  if  (B.9)  has  no  bounded  solution 
with  Re  Xg  >  0  then  (vc/wc)  is  linearly  stable  for  0  <  c  <<  1.  We  remark  that 
eigenvalue  problems  of  the  form  (B.9),  (B.10)  were  also  encountered  in  [16]. 

Equation  (B.9)  always  has  the  solution  XQ  =  0  ,  i.e.,  =  -1/y  /  and  vg = ^0  * 

This  corresponds  to  the  translation  invariance  of  ^Q( z)  and  does  not  determine 
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a  second  solution  to 


l  . 


!  ^ 

[ 

i 


stability.  Because  of  the  quadratic  relation  between  pQ  and  X^  , 

2 

{ B. 9 )  with  u0  *  -1/y1  ,  v0  -  *0  is  XQ  *  -eYj  +  1/y1  .  If  ty  <  1  then  this  second 

solution  satisfies  XQ  >  0  so  that  vc(z),  wc(z)  is  unstable.  On  the  other  hand,  suppose 
2 

ey  >1.  Then  the  second  solution  satisfies  X^  <  0  and  does  not  lead  to  an  unstable 
mode.  Moreover,  since  the  eigenfunction  v?q  is  of  one  sign  it  corresponds  to  the  lowest 
value  of  pg  in  the  discrete  spectrum  of  (B.9).  Therefore,  any  other  p^  must  be  greater 

than  — 1 / y^  and,  from  (B.10),  the  corresponding  values  of  Xq  satisfy  Re  XQ  <  0  .  Hence, 
2 

if  cy  >  1  we  have  linear  stability  for  small  c  • 


l 


! 


Without  loss  of  generality  we  take  xQ  -  0  .  If  yQ  <  ^  <  9  i,e*  1  <  2a/^v  <  2' 

then  for  PWL  the  steady  state  solution  to  (4*2)  is 

q(x)  -  a  ea*X~^  f  x  <  5 

(C.1) 

-  E  +  (a  -  E  /2)e0(*-e)  -  (E  .  e  <  *  <  0 


where 

(C. 2 )  C  -  TI  *n  (2a/E  -1) 

10  v 

2 

and,  as  in  Section  3,  a  m  1  +  1/y.  For  stability,  we  consider  exponential  solutions 
e*fc  v(x)  e*fcW(x)  to  the  linear  variational  equation*  As  in  Appendix  B,  W  satisfies 
W(x>  *  e  V(x)/(X  +  ey)  and  X,  V(x)  solve  the  Schroedinger  eigenvalue  problem: 

(C.3)  V"  +  [y  +  f * (q( x) )  ]  V  -  0  —  <  x  <  0 

where 

{ C. 4 )  y  *  y(  X)  *  -X  -  e/(  X  +  ey) 

and 

(C.5)  V(0)  -  0  ,  1  V(x)  |  V  0  as  x  X  . 

For  PWL, 

(C.6)  f  *  (q(  x) )  -  -1  +  6<5)/q'U>  • 

Here  6(5)  is  the  Dirac  delta  function. 

If  (C.3)  -  (C.6)  has  a  solution  with  Re  X  >  0  ,  then  we  say  that  q(x)  is 
unstable.  On  the  other  hand,  if  there  is  no  solution  with  Re  X  >  0  then  q(x)  is 
linearly  stable.  Indeed,  it  is  possible  that  there  may  be  no  values  of  X  for  which 
(C.3)- (C.6)  has  a  solution;  i.e.,  there  may  be  no  eigenvlaue  in  the  discrete  spectrum  for 
which  V(0)  *  0  .  In  such  a  case  we  would  conclude  that  q(x)  is  linearly  stable.  We 
note  in  passing  that  there  is  always  a  solution  y^,  Vg(x)  which  satisfies  (C.3)  and  for 
which  Vq(0)  *  0  ;  the  lowest  such  value  for  yQ  constitutes  the  "ground  state".  For  the 
solutions  we  seek,  the  relation  (C.4)  implies  that  the  largest  Re  X  corresponds  to  the 
algebraically  smallest  y  ;  for  this  smallest  y  ,  V  is  of  one  sign  for  x  <  0  . 


First  we  show/  for  "cubic- like-  f  ,  that  q(x)  is  linearly  stable  if  ey  >  1  •  If 
there  is  an  eigensolution,  it  satisfies  (C.3)  rewritten  as 

v-  +  [p  +  1/y  +  f '  (q( x) )  -  1/y}V  -  0  . 

Now  multiply  this  equation  by  q'(x)  and  (4*2)  by  V*(x),  add  the  two  equations,  and  then 
integrate  from  x  *  -•  to  x  *  0  to  obtain: 


q* (O)V'(O)  -  /°  ~  C-Vf(q)  +  Vq/y]  dx  -  (p  +  1/y)  f°  q'  V  dx  , 


Since  V(0)  *  0  and  f(Ev)  *  Ey/y  M  arrive  at 


M  +  1/Y  -  -q' (O)V'(O)/  f°  q'  V  dx  . 


But  q'fxJ/q'tO)  >  0  and  V(x)/V*(0)  <  0  for  x  <  0  .  Therefore  p  >  -1/y  ,  and  from 

2 

(C.4)  we  conclude  that  Re  X  <  0  when  cy  >  1  . 

2 

Next  we  consider  ey  <  1  and  restrict  attention  to  FWL.  If  we  normalize  V  so  that 
V(£)  -  1  ,  then 


V (x)  -  e 


a(x-£> 


*  sinh  opc/sinh  a£  ,  £  <  x  <  0 


Since  V  must  satisfy  the  jump  condition 


V’]5  -  -  V(5)/q'  u) 

C 

we  find  that  p  solves  the  "characteristic  equation w 
(C.8)  G(p,a,e,y)  «  0 

where 

( C .9 )  G(p,a,e,y)  =  1  -  coth  a£  -  1/(oao). 


Equation  (C.8)  can  be  analyzed  graphically.  For  all  values  of  the  parameters 
a,  e,y,  we  have  that  G  >  0  for  p  «  0  •  Next,  we  find  parameter  ranges  such  that  G  > 

0  for  0  <  1  -p  <<  1  •  In  such  cases  we  conclude  that  G  is  of  one  sign  for  p  <  1  • 

(If  not,  then  G  would  have  two  zeroes  for  p  <  1  .  This  would  imply  that  (C.3)-(C.6) 
has  t wo  different  eigenvalues  p  ,  each  with  V(x)  non-zero  for  x  <  0  ,  which 
contradicts  the  known  spectral  properties  of  Schroedinger  problems.)  From  this  it  follows 
that  there  is  no  solution  to  (C.3)-(C.6)  of  the  assumed  form  so  that  q(x)  is  linearly 
stable  for  such  parameter  ranges. 

To  find  the  parameter  range,  suppose  p  is  slightly  less  than  one.  Then  G  >  0  is 

equivalent  to  -£  <  aa  or  ,  using  (C.2) , 

(C.10)  -in(2a/E  -  1)  <  2a/E  . 

v  v 

2 

Hence,  if  ey  <  1  and 

2a/Ev  >  1.278...  , 

the  inequality  (C.10)  is  satisfied,  there  is  no  eigensolution  of  the  desired  form,  and 
therefore,  q(x)  is  linearly  stable. 

To  see  when  q(x)  is  unstable,  observe  from  (C.4)  that  p  <  -ey  implies  that  at 
least  one  of  the  two  corresponding  X  values  has  Re  X  >  0  •  Hence,  if  G  <  0  for 
p  =  -ey  ,  then  G  has  a  zero  for  p  <  -ey  and  q(x)  is  unstable.  For  fixed 
e,y  and  p  *  -ey  ,  we  write  G  <  0  as 
(C.11)  1  -  coth [ / 1  +  ey  Jin  (2a/Ev  -  1)/2d  <  2a/[/l  +  ey  2a/EvJ . 

If  0  <  2a /E  -  1  <<  1  ,  (C.11)  becomes  1  <  o//l  +  ey  which  is  satisfied  when  ey^  <  1 
so  that  q(x)  is  unstable.  Finally,  observe  that  at  critical  values  of  the  parameters 

a,e,y  for  which  G(-ey,a,e,y)  ■  0  ,  the  corresponding  X  values  are  a  complex  pair  with 
j  2~ 

Re  X  *  0,  la  X  *  i  i  /e(1  -  ey  )  so  that  a  branch  of  periodic  solutions  emerges  through 
Hopf  bifurcation. 
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